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ABSTRACT 

Photons  are  seriously  scattered  when  entering  turbid  medium;  thus  the  images  of  objects  hidden  in  turbid  medium  can  not  be 
obtained  by  just  collecting  the  transmitted  photons.  Early-arriving  photons,  which  are  also  called  ballistic  or  snake  photons, 
are  much  less  scattered  when  passing  through  turbid  medium,  and  contains  more  image  information  than  the  late-arriving 
ones.  Therefore,  objects  embedded  in  turbid  medium  can  be  imaged  by  gathering  the  ballistic  and  snake  photons.  In  the 
present  research,  we  try  to  recover  images  of  objects  in  turbid  medium  by  simultaneously  “time-gate”  and  “polarization- 
gate”  to  obtain  the  snake  photons.  An  Argon-pumped  Ti-Sapphire  laser  with  lOOfs  pulses  was  employed  as  light  source.  A 
streak  camera  with  2ps  temporal  resolution  was  used  to  extract  the  ballistic  and  snake  photons.  Two  pieces  of  lean  swine 
meat,  measured  4mmX3mm  and  5mmX4mm,  respectively,  were  placed  in  a 10cmX10cmX3cm  acrylic  tank,  which  was  full 
of  diluted  milk.  A pair  of  a polarizer  and  an  analyzer  was  used  to  extract  the  light  that  keeps  polarization  unchanged.  The 
combination  of  time  gating  and  polarization  gating  resulted  in  good  images  of  objects  hidden  in  turbid  medium. 


1.  INTRODUCTION 

1.1  Recent  development  of  near-infrared  light  imaging 

Recently,  optical  imaging  through  highly  scattering  medium  with  near-infrared  light  has  received  much  attention  [1-8]. 
Because  of  its  noninvasive  nature,  researchers  have  regarded  near-infrared  light  as  a safe  alternative  to  ionizing  radiation 
such  as  x-ray  or  positron  emission  tomography.  However,  near-infrared  light  can  penetrate  into  biological  tissues  only  a 
limited  depth  because  of  strong  scatter.  Therefore,  advanced  optical  techniques  are  required  for  such  imaging  application. 

It  was  shown  that  near-infrared  light  is  only  weakly  absorbed  by  most  biological  tissues,  and  light  propagating  in  biological 
tissues  is  mainly  attenuated  through  scattering  [9-11].  Based  on  this  fact,  two  approaches  were  commonly  adopted  for 
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imaging  with  near-infrared  light:  1)  capturing  snake  photons  for  direct  imaging  [2,8,12,13],  and  2)  collecting  diffusive 
photons  and  reconstructing  the  image  with  the  theory  of  diffuse  photon  density  wave  [1,5,10,11].  The  second  approach  costs 
less  than  that  of  the  first  one;  however,  its  spatial  resolution  is  usually  unsatisfied.  Due  to  multiple  scattering,  the  image 
information  of  the  embedded  objects  for  the  diffusive  light  are  almost  lost.  In  fact,  the  spatial  resolution  of  reconstructed 
images  by  using  diffusive  photons  is  limited  by  many  factors,  including  optical  properties  of  the  turbid  medium,  the 
thickness  of  the  samples,  the  wavelength  of  the  probe  light,  etc.  Therefore,  the  reconstruction  of  high-resolution  images  is 
not  a cakewalk  [14].  On  the  other  hand,  the  collection  of  ballistic  or  snake  photons  can  provide  high-resolution  image,  but 
the  amount  of  these  photons  are  very  small  and  not  easily  captured.  Furthermore,  these  photons  exists  in  a very  short  time 
period  (200ps~300ps),  dependent  on  the  optical  property  of  the  measured  tissues.  In  recent  years,  due  to  the  development  of 
high  speed  detector,  such  as  electronic  gated  CCD  and  streak  camera,  the  research  on  imaging  techniques  with  ultrashort 
pulses  has  gradually  become  popular  [13,15]. 

1.2  polarization  gating  and  time  gating  for  snake  photons 

Since  snake  photons  encounter  little  scattering  and  are  along  shorter  path  when  transmitting  the  measured  sample,  they 
arrive  at  the  detector  earlier  than  diffuse  photons.  By  gating  those  early  arriving  photons,  we  are  able  to  obtain  the  images  of 
the  embedded  objects.  However,  usually  it  is  not  easy  to  determine  the  time  interval  of  snake  photons.  Another  property  of 
the  snake  photons  is  that  their  polarization  almost  keeps  unchanged  during  propagation.  In  other  words,  if  the  incident  pulse 
is  polarized,  the  output  snake  photons  bear  the  same  polarization  [16].  Therefore,  combining  time  gating  and  polarization 
gating  is  useful  for  extracting  snake-photons  and  get  clear  images. 


2.  EXPERIMENTAL  PROCEDURES 

2.1  Samples — turbid  medium  and  targets 

In  the  experiment,  diluted  milk  (the  milk  was  made  up  by  milk  power  branded  “Kolin”),  was  used  to  simulate  turbid 
medium.  The  milk  was  contained  in  an  transparent  acrylic  tank,  with  10cmX10cmX3cm  in  dimensions.  Two  pieces  of  lean 
swine  meat,  measured  4 mm><3  mm  and  5 mm><4  mm,  respectively,  hanged  in  the  center  of  the  tank  and  separated  by  about  3 
mm,  were  embedded  in  the  milk  as  objects  to  be  imaged.  The  thickness  of  the  both  meat  was  about  1 mm.  The  total 
attenuation  coefficient  (//,)  of  the  objects  and  the  background  medium  are  13.6cw-1  and  10 cm~] , respectively.  The  tank 
was  then  placed  on  a motor-controlled  two-dimensional  translation  stage  for  measurement. 

2.2  Measurement  instruments 

The  experimental  setup  is  shown  in  Fig.  1.  An  Argon-pumped  Ti-Sapphire  laser  (Spectra  Physics,  USA),  which  emits  lOOfs 
short  pulses,  was  employed  as  the  light  source.  The  streak  camera  was  operated  in  the  synchronous  scan  mode,  which  would 
accumulate  the  light  intensity  of  many  pulses,  to  obtain  a high  signal  to  noise  ratio.  One  of  the  split  light  from  the  first  beam 
splitter  (BSi)  was  incident  on  a PIN  diode  to  generate  a trigger  signal.  The  other  part  of  light  was  separated  into  two  paths 
by  the  second  beam  splitter  (BS2).  One  is  then  incident  on  the  samples  for  imaging,  while  the  other  was  sent  to  the  streak 
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Fig.l  Experimental  setup. 


camera  (Hamamatsu,  K.K.,  Japan)  directly  as  the  reference  signal.  During  measurement,  the  incident  light  and  the  detector 
were  fixed  and  aligned  with  each  other,  while  the  tank  was  moved  with  a two-dimensional  motor-controlled  translation  stage. 
The  measurement  area  was  2cmX2cm  in  dimensions.  The  temporal  profiles  of  pulses,  which  were  broadened  (diffused)  to 
become  about  1.1  ns  in  width  after  transmitting  the  turbid  medium,  were  collected  point  by  point.  Two  polarizers  Pj  and  P2, 
placed  before  and  after  sample,  were  used  to  gating  the  output  photons  whose  polarization  remained  unchanged.  This  will 
help  the  streak  camera  extracting  the  early-arriving  photons,  whose  polarization  was  supposed  to  be  unchanged. 

2.3  Instrument  control  and  data  acquisition 

The  instrument  control  and  data  acquisition  were  implemented  by  using  software  named  LabVIEW  (National  Instruments 
Inc.,  USA).  The  movement  of  the  two-dimensional  translation  stage,  including  the  moving  speed,  step,  displacement,  axis, 
etc,  can  be  controlled  by  the  self-developed  software.  The  timing  for  sample  moving  was  pre-determined  to  match  time 
required  for  getting  profiles  by  the  streak  camera.  The  acquired  temporal  profile  was  stored  and  translated  into  a text-format 
file.  The  data  processing,  including  intensity  normalization,  temporal  calibration,  etc,  was  completed  with  an  analysis 
program  written  in  LabVIEW.  After  calculation,  the  data  were  stored  in  the  hard  disk,  and  the  image  was  displayed  on  the 
screen. 
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3.  IMAGING  RESULTS 


3.1  Temporal  profile 

Temporal  profiles  of  output  pulses  can  be  obtained  by  integrating  the  streak  image  for  each  measuring  point  of  the  sample. 
The  reference  beam,  which  did  not  pass  through  the  turbid  medium  and  always  arrived  at  the  detector  periodically  with 
fixed  intensity,  was  used  for  time  and  intensity  calibration.  Figure  2 shows  a typical  temporal  profile.  Occasionally  jitters 
occur  due  to  the  instability  of  the  pulse  repetition  rate  and  laser  intensity.  To  remove  such  undesired  effects  we  normalize  the 
intensity  of  the  signals  based  on  that  of  the  reference.  Also,  the  temporal  positions  of  the  profile  were  adjusted  according  to 
the  center  of  reference.  Figure  3 shows  a typical  calibrated  temporal  profile.  In  this  chart,  all  data  points  were  shifted  such 
that  the  peak  of  the  reference  signal  was  positioned  at  zero  time  (t  = 0). 

3.2  Measurement  with  polarizers 

To  further  reject  the  diffusive  photons  for  improving  the  imaging  results,  two  polarizers  were  placed  before  and  behind  the 
sample  (see  Fig.  1).  The  first  polarizer  (PJ,  which  was  placed  before  the  sample,  was  used  to  generate  a linearly  polarized 
incident  light.  The  second  one  (P2),  which  was  placed  behind  the  sample,  was  used  as  a polarization  analyzer.  During  the 
measurement,  the  analyzer  (P2)  was  first  rotated  to  have  the  same  orientation  of  P„  and  the  measured  signal  was  denoted  as 
Lp  (parallel  polarization).  Next,  it  was  rotated  to  become  perpendicular  to  P„  and  the  measured  signal  was  denoted  as  Lc 
(cross  polarization).  It  is  supposed  that  ballistic  and  snake  photons,  which  are  almost  not  scattered  while  propagating  in  the 
turbid  medium,  maintain  the  original  polarization.  However,  the  diffusive  photons,  which  has  been  seriously  scattered,  have 
lost  the  original  polarization  information  after  transmitting  through  the  scattering  medium,  and  the  orientations  of 
polarization  for  such  photons  were  random.  About  50%  of  the  diffusive  light  passes  the  parallel-  and  cross-polarization 
analyzer,  respectively.  Hence,  the  summation  of  Lp  and  Lc  represents  the  total  outptu  light.  Only  those  photons  keeping 
their  original  polarization  left  in  the  profile  of  Lp-Lc.  Figure  4 shows  the  results  for  summation  and  subtraction  of  Lp  and  Lc. 
In  this  figure,  we  saw  that  the  diffusive  light  (the  late-arriving  photons)  still  exists  in  Lp+Lc,  while  only  the  ballistic  and 
snake  photons  (close  to  the  leading  edge  of  the  curve)  left  in  the  curve  Lp-Lc.  Therefore,  we  conclude  that  with  two 
polarizers  and  a streak  camera,  we  can  extract  the  light  keeping  the  original  polarization,  which  bears  abundant  information 
for  imaging. 

To  evaluate  the  possibility  for  optical  imaging  with  the  described  method,  we  scanned  a 2cmX2cm  area  for  images.  The 
depth  of  the  tank  was  about  3cm.  For  each  measuring  point,  the  parallel  polarized  light  (Lp)  and  cross  polarized  light  (Lc) 
were  both  measured.  Figure  5 shows  the  image  that  was  obtained  from  Lp+Lc,  the  total  light  passing  the  turbid  medium. 
Since  the  signal  contains  diffusive  photons,  we  can  hardly  recognize  the  shape  of  the  targets. 

On  the  other  hand,  Fig.  6 shows  the  image  of  obtained  from  the  profile  Lp-Lc.  As  the  inset  shows,  only  the  early  arriving 
photons  (near  the  leading  edge  of  the  profile)  are  used.  Further  calculation  shows  that  the  contract  (defined  as  the  image 
minus  background  to  the  background  ratio)  of  the  image  in  Fig.  6 are  about  2.5  times  of  that  in  Fig.  5.  In  other  words,  in  Fig. 
6 the  embedded  objects  can  be  identified  more  clearly  than  in  Fig.  5.  The  white  square  boxes  indicate  the  actual  positions  of 
the  two  pieces  of  meat. 
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4.  DISCUSSION  AND  CONCLUSION 


In  Fig.  6,  the  left  image  is  much  clearer  than  the  right  one.  In  fact,  the  right  one  contains  more  fat  than  the  left  one.  Since  the 
color  of  fat  is  almost  white  at  room  temperature,  the  scattering  effect  of  this  piece  of  meat  was  much  severe  than  regular  lean 
meat.  Therefore,  the  image  of  the  right  one  became  blurrier  than  the  left  one. 

In  this  research,  by  combining  time  gating  and  polarization  gating,  we  obtained  the  images  of  biological  tissues  embedded  in 
turbid  medium  (diluted  milk).  It  was  shown  that  polarization  gating  was  an  efficient  approach  that  could  help  extracting 
ballistic  and  snake  photons.  Therefore,  clearer  image  can  be  obtained.  Furthermore,  the  addition  of  polarization  gating 
provides  a useful  criterion  forjudging  the  interval  of  “snake  photons”,  which  would  be  very  useful  for  the  study  of  optical 
properties  of  tissues  to  be  measured. 


Fig.  4 Temporal  profiles  for  the  summation  of  parallel-  and  cross-polarization  light  (Lp+Lc,  which  represents  the  photons 
passing  through  the  sample  and  analyzer  without  considering  the  polarization),  and  the  subtraction  of  parallel-  and 
cross-polarization  light  (Lp-Lc,  which  represents  the  photons  without  changing  their  polarization  after  passing 
through  the  sample). 
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